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ABSTRACT 

QSOs have been thought to be important for tracing highly biased regions in the early universe, 
from which the present-day massive galaxies and galaxy clusters formed. While overdensities of star- 
forming galaxies have been found around QSOs at 2 < z < 5, the case for excess galaxy clustering 
around QSOs at z > 6 is less clear. Previous studies with HST have reported the detection of small 
excesses of faint dropout galaxies in some QSO fields, but these surveys probed a relatively small 
region surrounding the QSOs. To overcome this problem, we have observed the most distant QSO at 
z = 6.4 using the large field of view of the Suprime-Cam (34' x 27'). Newly-installed red-sensitive fully 
depleted CCDs allowed us to select Lyman break galaxies (LBG) at z ~ 6.4 more efficiently. We found 
seven LBGs in the QSO field, whereas only one exists in a comparison field. The significance of this 
apparent excess is difficult to quantify without spectroscopic confirmation and additional control fields. 
The Poisson probability to find seven objects when one expects four is ~10%, while the probability 
to find seven objects in one field and only one in the other is less than 0.4%, suggesting that the 
QSO field is significantly overdense relative to the control field. These conclusions are supported by a 
comparison with a cosmological SPH simulation which includes the higher order clustering of galaxies. 
We find some evidence that the LBGs are distributed in a ring-like shape centered on the QSO with 
a radius of ^3 Mpc. There are no candidate LBGs within 2 Mpc from the QSO, i.e., galaxies are 
clustered around the QSO but appear to avoid the very center. These results suggest that the QSO 
is embedded in an overdense region when defined on a sufficiently large scale (i.e. larger than an 
HST/ACS pointing). This suggests that the QSO was indeed born in a massive halo. The central 
deficit of galaxies may indicate that (1) the strong UV radiation from the QSO suppressed galaxy 
formation in its vicinity, or (2) that star-formation closest to the QSO occurs mostly in an obscured 
mode that is missed by our UV selection. 



1. INTRODUCTION 

At 2 < z < 5, strong overdensities of star-forming 
galaxies have been found around QSOs and radio galax- 
ies, and thus, QSOs have been thought to trace highly 
biased regions in wh ich the present-day massive galaxy 
clusters formed (e.g. iDiorgovski et"aTJl2003t iMilev et"alj 
120041 : [Kashika wa et al.l l2007). It is expected that this ex- 
tends to the most luminous z ~ 6 QSOs as well, as this 
rare population hosts supermassive black holes of several 
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billion solar masses that are presumed to reside in the 
most massive galaxies and dark matter halos present at 
this redshift. 

However observational results to date have been puz- 
zling; Five z ~ 6 QSO fields observed by the HST/ACS 
show no major enhancements in the galaxy density 
(|Kim et al.ll2009t ). Even though a few QSO field s have 
showed an appar ent overdensity (|Stiavelli et al.l 120051 : 
iZheng et al.1 l2006h . none of them were among the rich- 
est structures discovered to date at z ~ 6, as evi- 
denced by much larger overdensities found in random 
fields (e.g., lOuchi et alj|2004bt iKashikawa et al.l [20071 
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lOta et all [20081 see Overzier et al. (2009) for a discus- 
sion) . 

Why do QSOs suddenly appear to stop being at the 
center of the overdensity at z ~ 6? One hypoth- 
esis is that despite the higher dark matter densities 
near the QSOs, the strong ionizing UV radiation from 
the central QSOs may prohibit the condensation of gas 
thereby suppressing gal axy formation around the QSOs 
(|Barkana fc Loebl [1999) . An alternative and perhaps 
more likely explanation is that the lack of overdensities 
identified is related to the fact that it is currently tech- 
nically challenging to perform a survey deep enough to 
detect faint z ~ 6 galaxies and cover an area large enough 
to probe the large-scale structure surrounding the QSOs. 

In this work, we aim to study the large scale structure 
around the currently most distant QSO at z = 6.43, tak- 
ing advantage of the large aperture of the 8.2m telescope 
"Subaru" located on Maunakea Observatory and the 
wide field prime focus c amera "Suprime-Cam" (34'x27') 
(|Mivazaki et al.ll200 2b). In addition, we have recently 
installed new red-sensitive fully d epleted CCDs on the 
Suprime-Cam (jKamata et al.ll2008l ) , allowing us to reach 
necessary depths in much shorter exposure time. Im- 
proved sensitivity is a factor of about 1.4 and 2 better in 
the z'-band and at ~ l//m, respectively. 

Throughout the paper, we use Hq = 70/i7o km s _1 
Mpc -1 , D, m = 0.3, ft A = 0.7. Magnitudes are given m 
the AB system. 

2. DATA & ANALYSIS 

2.1. Observation 

We observed the field surrounding the most distant 
QSO, CFHQS J2329-0301 at z = 6.43, M U5 o = -25.23 
(Will ott et al.ll2007h using Suprime-Cam with the filters 
i',z', zr during an engineering run in Aug 2008 and UH 
time in Jun 2 009. We used the speci al zr band filter 
constructed by iShimasaku et al.l (|2005l ) . This filter cov- 
ers the redder side of the z' band and has a central wave- 
length of 9900A (Figfrj . The seeing was stable through- 
out the runs, with a FWHM w 0.5 arcsec in i' , 0.4 ~ 0.7 
arcsec in z' and 0.5 ~ 0.7 arcsec in zr. The exposure 
times were 360 and 500 s in i' and zr, respectively. The 
exposure times in the z' band varied; we obtained 3 ex- 
posures of 100 s, 2 of 300 s, 3 of 400 s, 6 of 500, and 2 of 
700 s resulting in more than 10000 ADU of sky counts. 
The total exposure time in i', z' and zr were 3,600 s on 
Aug 28 2008, 6,900 s on 2, Aug 27 and 12,532 s on Aug 
27 2008 and Jun 18 2009 respectively, under good see- 
ing conditions (less than 0.6 arcsec FWHM). Our dither 
strategy consisted of five or more pointings on a circle of 
1 arcmin radius. 

2.2. Data Reduction 

Our reduction procedure follows iMivazaki et al.l 
(2002a), with small adjustments made related to the 
new CCD. A major difference between the old and new 
Suprime-Cam is that the latter has 4 channel read out 
circuits for faster readout. As a result, over-scan sub- 
traction and background sky-subtraction should be per- 
formed on each channel individually for each CCD. An 
overscan value is evaluated by taking the mean along the 
columns of the overscan region of each channel, and then 
subtracted from the science region. A flat field image is 



created by taking the median of all science images after 
masking objects. The background sky is subtracted on 
each channel after flat fielding. To estimate the back- 
ground, we first reject data points having 2 times larger 
or 4 times lower sigma value than the sky value. The 
upper cut prevents astronomical objects from contam- 
inating the sky, while the lower cut removes residuals 
from the flat-fielding. Then, we estimate the modes in 
64 pixel x 64 pixel boxes. A background of each pixel is 
estimated by interpolating these mode values at vertexes 
of a triangle centered on the pixel. 

Before combining all the science frames, the images 
need to be corrected for (1) geometric distortion, (2) 
displacement and rotation of each CCD from the fidu- 
cial position, and (3) pointing offsets between exposures. 
We corrected for these effects by minimizing the posi- 
tional differences of common control stars identified on 
each CCD in all exposures relative to the first exposure. 
The resulting positional error of the control stars from 
the fiducial point is approximately 0.5 pixel (rms). We 
solve the transformation among the 3 bands, aligning 
each frame with the first exposure of the zr band in or- 
der to perform multi-band photometry using an aperture 
as small as possible to obtain a high S/N ratio. Flux off- 
sets are also calculated during this procedure. Since we 
want to determine the transformation among the bands 
as accurately as possible, an image warping procedure 
was performed only once, and a second order bi-linear 
interpolation was used for rebinning. For stacking, we 
adopted the clipped mean in order to eliminate cosmic 
rays. The implementation of these reduction steps were 
performed using the software suite imcat. 

2.3. Photometric Calibration 

Photometric calibration is performed by com- 
paring stars in a reference field 2 deg away 
to the North (a,6) 20 oo = (23 : 29,-01 : 01) 
with those in the SDSS/DR7 s tar catalog 
(jAbazaiian fc Sloan Digital Sky Survey! 120081 ) and 
using the following equations: 

i'subam = 0.125(i' - r') SDSS + 0.003 + i' SDSS (1) 
z'subam = -1.091(i' - z') SDSS + 0.004 + i' SDSS (2) 
ZR,Subaru = -1.414(i' - z') S dss + 0.021 + i' SDSS (3) 

These equations are determined from iGunn fc Strykerl 
(1983) and convolved with response curves that include 
both optics and the atmos phere. We mea s ured the effi- 
ciency of the new CCDs in iKamata etaLI (|200l . since 
the calibration field was observed soon after the com- 
pletion of our science exposure, no airmass/ atmospheric 
corrections were needed to obtain the photometric zero- 
point. After the correction, the rms of the magnitude 
differences between our catalog and the SDSS catalog 
are 0.06, 0.07 and 0.07 in z' and zr, respectively. 

We have checked the colors of a sample of star-like 
objects selected according to CL ASS_STAR>0 . 9, by com - 
paring the color of the stars from lGunn fc Strykerl (|1 9831 ) 
to the observed ones, confirming a good internal consis- 
tency of colors in our catalog. The template and ob- 
served colors were consistent within an accuracy of 0.05 
mag. Since the seeing in the final images are is 0.58, 0.54 
and 0.50 arcsec (FWHM), we use a small aperture to 
perform photometry to obtain a better S/N ratio. Fol- 
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lowing iShimasaku et al.l (|2005l ). who used an aperture 
size twice that of PSF, we adopt a 1.2 arcsec aperture to 
perform photometry. The 3 sigma limiting magnitudes 
within this 1.2 arcsec aperture are i' = 26.95, z' = 26.13, 
and zr — 25.46 (AB) mag. Here, 1 a sky magnitudes 
are computed by randomly placing a 1. 2 arcsec aperture 
in the blank (sky) part of each image (lYaei et all 1200a 
IQuchi et al.ll2004al) . 

2.4. Astrometric Calibration 

We obtain an astrometric solution by compari ng our 
catal og with the USNO-B1.0 catalog (see scamp (jBertinl 
2006)). The resulting accuracy of the astrometric solu- 
tion is 0.5 arcsec (rms). Note that we did not use this 
astrometric solution to stack individual raw images, but 
only to obtain the final absolute positions of the objects. 

2.5. Object Detection 

The object detection is performed by SExtractor 2.3.2 
(jBertin fc Arnouti I1996D . We constructed a detection 
image, the "&W-ZR image, by combining all science 
frames, including those with a slightly poorer seeing con- 
ditions. The object detection reliability using this "all"- 
zr image is higher than when only using the zr image 
for detection given that w 3cr sources in the zr image 
are detected at w 3.7a in the "&[1"-zr image. Using the 
"all"-Zfl image, we consider an object detected if it has 
more than 5 connected pixels with each exceeding the lo- 
cal sky rms by a factor of 2. To reduce contamination by 
false detections, we only use the "detected cleanly" and 
"no blending" objects that have FLAGS=0 in SExtractor. 
We measure the magnitudes in 1.2 arcsec apertures for 
each pass-band to derive colors of the detected objects, 
and adopted MAG_AUT0 as our estimate for the total zr 
magnitude. Objects in the ?'-band fainter than 2 sigma 
were replaced by the 2 sigma limiting magnitude as an 
upper limit. In addition, all objects having detections 
in the z'-band of less than 2 a were rejected in order to 
keep spurious detection at a minimum. 

Some lower S/N regions, such as the "blooming re- 
gions" (the halos and horizontal spikes surrounding 
bright stars) , as well as the outer edge of the image, are 
masked manually. The resulting effective field of view is 
0.219 deg 2 . We detect 48,632 objects after masking (to a 
limiting magnitude of 25.46, or 3.7 sigma detections on 
the "all"-ZR image). The number of detections is com- 
parable to that in our comparison field, the Subaru deep 
field (SDF), which contains 45 ,405 objects in a single field 
of view of the Suprime-Cam ([Shimasaku et a l. 2005). 

2.6. SED Modeling 

We have upgraded the CCDs of the Suprime-Cam to 
red-sensitive ones, which have 1.4 and 2 times better sen- 
sitivity in z' and zr compared to the previous detectors. 
Thus, we can reach the necessary depth in much shorter 
exposure time. This improved sensitivity allows us to 
efficiently use the z' — zr color to select Lyman break 
galaxies at z > 6.4. At the QSO redshift of z = 6.43, the 
z' — zr color selection is more efficient in selecting galax- 
ies than the i' — z' color selection (the so-called i-dropout 
technique) , which has been used in previous works to se- 
lect 5.5 < z < 6.5 galaxies (FigfTJ). 

In order to predict the colors of galaxies as a function of 
their redshift, we computed z' — zr colors of model galax- 
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Fig. 1. — An example of the spectra of a galaxy and a QSO 
redshifted to z = 6.43. The filter response curves in i' . z' and 
zr are indicated from left to right. The spectrum of the galaxy is 
based on a model having an exponentially-declini ng star formation 
history and an age of 0.5 Gyr (Bruzual & Chariot 2003), while the 
spectrum o f the QSO is based on the QSO composite spectrum 
taken from l|Vandcn Berk . D. E. et alJl20CT . 
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Fig. 2. — i' — z' vs z' — zr color-color diagram, showing the 
locations of detected objects, candidate z ~ 6.4 galaxies, stellar 
objects, and tracks indicating the expected colors of low- and high 
redshift galaxies. Gray dots are all objects fainter than zr = 24.0, 
Red squares are z ~ 6.4 galaxies selected in the QSO field. Blue 
plusse s are objects detected as z ss 6 galaxies in Shimasaku ct al. 
(2005) in the SDF. Gray box crosses are stars IIGunn &; S trvkcr 
ll98al(Jhiu et al.l2006HGolimowski et al.l20"0llKnapp et al.l2004f r 
The purple, green, and orange lines represent model galaxies hav- 
ing star "instantane ous burst", "exponential" and "constant" star 
formation histories (Bruzual & Chariot 2003). Absorption by the 
neutral hydroge n at each red shift was added following the pre- 
scription of Fan ct al. (2006a). The blue arrow indicates CFHQS 
J2329-0301 (see Table [Jk 



ies (0 .05, 0.1, 0.2 Gyr of age, using [Bruzual fc Charlotl 
(2003)) and added the absorption effec ts from the neutra l 
hydrogen in the intergalactic medium (jFan et al. 2006a) . 
According to our modeling, colors of i — z' > 1.9 and 
z' — zr > 0.3 can be used to identify galaxies at z > 5.8 
(Fig. 0. 

2.7. Object Selection 

Here we will select LBGs at z > 6.4 using our i', z', zr 
catalog. The model colors shown in Fig. [5] predicts that 
galaxies at z = 6.4 should have a large color difference of 
about i'—z' ~ 3.5. However, for an object having z' ~ 25 
this would require i' ~ 28.5 mag, much fainter than our 
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detection limits. Therefore, we slightly loosen the color 
cut and require %' — z' > 1.9 and z' — zr > 0.3. These 
relaxed criteria are still appropriate for selecting LBGs 
at z > 6.4, although we must be somewhat cautious of 
potential galaxies at z ~ 1.8 having colors ofi' — z' ~ 1.9 
that could make it into the sample as well. In order to 
assess the impact of this effect on our conclusions, in 
section [331 below we will try to statistically estimate the 
number of such z ~ 1.8 interlopers expected. 

In our control field (the SDF), only a single object was 
found that satisfied our color criteria, a nd this object was 
classifi ed as a genuine z ~ 6 galaxy by iShimasaku et ail 
(2005). Because the SDF data are deeper than ours, we 
adjusted the depth of that field by replacing the mag- 
nitudes of all objects fainter than our limits with the 
limiting magnitudes as computed for our field. We also 
removed two objects brighter than z^Auto = 24.0 since 
such objects are too bright to be at z = 6.4 (according to 
az~6 LBG study in the SD F there is no z > 6 galax y 
brighter than Zr — 24.8, see IShimasaku et "all (|20Q5T> ). 
These two objects may be contaminating stars since they 
have SExtractor CLASS_STAR values greater than 0.95, 
meaning that they are likely to be unresolved. 

As a result, to the magnitude limit of z' < 25.46, we 
have detected 7 objects that are good candidates for be- 
ing z ~ 6 LBGs. These objects are shown as red squares 
in Fig. [2j Note that the (lower limits on) their i — z 
colors are bluer than expected from the models, due to 
the relatively shallow limiting magnitude in i' . 

3. SOURCES OF CONTAMINATION 

3.1. False Detections 

Because we are selecting very faint objects, our sample 
may be contaminated by false detections due to back- 
ground fluctuations. To evaluate the number of such 
contaminations, we created a negative of the "all"-ZR im- 
age, and repeated our detection, photometry and mask- 
ing routines on each of the filter images. No detections 
were obtained based on this negative image. We conclude 
that our catalog is not affected by spurious detection. 

3.2. Contamination by Faint Dwarf Stars 

As we can see in Fig. [51 L/T dwarf stars also satisfy 
our adopted color cut (i' — z' > 1.9 & z' — zr > 0.3). In 
this subsection, we estimate the expected number of L 
and T dwarfs, given the size, depth and Galactic position 
of the field. 

As a first test, when we applied our color criteria to 
the SDF only one object was found and classified as a 
genuine z ~ 6 object (see above). The second reddest 
object in the field has i — z ~ 1.65, i.e., there are no 
stellar-like objects in the SDF with colors red enough to 
make it into our selection. 

However, we need to take into account the different 
galactic latitudes of our field and the SDF. Unfortu- 
nately, the late-type star counts as a function of Galac- 
tic position are not accurately known at the faint mag- 
nitudes we are probing. Late-type stars having similar 
colors to high-z galaxies can only be distinguished by 
using deep spectroscopic observations. The lat est es- 
timate was performed by iCaballero et al.l (|2008) . who 
studied the contamination from L, T dwarfs by mod- 
eling the spatial distribution of late type stars in the 



Galactic thin disc described using an exponential law 

m = n A<i e-^fm,^ = cosb^£sj_ ± ^ where nAi 

is the number density of i-type stars, Hr, hz are the 
scale length and the scale height of the exponential 
disc and d is the distance to us in galactic coordinates, 
d= (R% +ePcos 2 6- 2i? dcos6cosO 1/2 - The adopted 
parameters are R = 8.6kpc, Iir — 2.25kpc, hz — 0.3kpc. 
Local spatial densities of late type star s are a dopted from 
a compilation given inlCaballcro et al.l (|2008|) . while rela- 
tions bjtweeji^h^ypeand absolute magnitude are taken 
from iKnapp et all ([2004D . According to this model cal- 
culation, we expect 1.5 times as many L/T dwarfs as 
in the SDF for zr < 25.25 mag, corresponding to the 
faintest bin of the magnitude of our selected objects. 

We have seen that there were no stars in our z ~ 6.4 
color-selection in the SDF. Here we try to scale that null 
detection in the SDF to our field. If we assume that 
stars follow a Poisson distribution, the expected number 
of L/T dwarfs in the SDF is 0.40 per field at a 68% confi- 
dence level. Using the Galactic scaling computed above, 
we expect that the probability of finding two or less stars 
in our field is 87%. Given that we have already rejected 
two stellar like objects (§2.5), the chance of finding one 
more star is less than two%. Therefore, we conclude that 
there is little chance that our z ~ 6.4 LBG catalog is 
contaminated by stellar objects. We note, however, that 
object with id 1 could be a dwarf star since its i' — z' 
color is not large enough. We will include this object in 
the following discussion. 

3.3. Contamination by red galaxies at z ~ 1.8 

As mentioned earlier, we found only one object when 
applying our adopted color criteria for z^6.4 galaxies to 
the SDF catalog (after adjusting the limiting magnitudes 
to ours), and this object was identified as a z ~ 6 object 
given that its B, V, R magnitudes were too faint for a 
low-z interloper. Based on this, the expected number of 
z ~ 1.8 galaxies found according to our criteria should be 
close to zero with an upper limit of 1.83 at a confidence 
level of 84% (Gehrels 1986). Although cosmic variance 
should ideally be taken into account as well, we do not 
have any other suitable comparison fields that are as deep 
in zr as the SDF. 

How ever, we use the pub licly available SXDS DR1 cat- 
alog (jFurusawa et al.ll2008D to obtain a rough estimate of 
the expected contamination using only ai-z > 1.9 color 
cut. The SXDS covers 5 pointings with the Suprime- 
cam, and thus, the cosmic variance is less of a concern. 
We first applied again the observational limits as given 
by our field by replacing all magnitudes that are fainter 
than our limiting magnitudes (i' 2a . — 27.39, z' 2a . — 26.57) 
by our 2a upper limits. We used FLAG=0 in order to se- 
lect only cleanly detected objects, and trimmed a region 
of lower S /N (500 pixels from the edge of the field) . Only 
3 objects were found that passed the selection criteria in 
all 5 pointings of the SXDF. These objects are detected 
in the B band (brighter than 28.4, 3cr mag), indicating 
that they are most likely low redshift interlopers at a 
redshift of z ~ 1.8. If we scale this number to the area 
of our single pointing field, only 0.6 of such interlopers 
are galaxies are expected in our. This rough estimate is 
consistent with the estimate derived above based on the 
SDF, suggesting that the contamination from z ~ 1.8 
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Fig. 3. — Number counts of all detected objects and the z-dropout 
galaxies. We did not apply completeness corrections. The number 
of LBGs in our field is larger than that in the SDF although the 
opposite is true when comparing the total number of 1 detected 
objects. 



galaxies is negligible. 



4. RESULTS 



4.1. Number Counts 

Fig. |3] shows the number counts of all detected objects, 
as well as those of the z ~ 6 galaxies. The number of all 
objects in this field is slightly lower, by about 20%, than 
those detected in the SDF field at zr < 24. This could 
be due to a small difference in the absolute photometric 
calibrations between the two catalogs, or due to cosmic 
variance. If it were due to a difference in photometric 
calibration, a shift of 0.2 magnitude would be required 
to explain the difference. Note that such an offset, if real, 
would not affect the colors used in our color selection of 
the z ~ 6 galaxies, as the colors were checked for inter- 
nal consistency using a catalog of stars . We next con- 
sider t he possibility of cosmic variance. iFurusawa et al.l 
( 2008) observed five Suprime-Cam pointings and derived 
the number counts. The variance in the number counts 
was found to amount to a factor of 1.7 among the in- 
dividual pointings in the z' band. Thus, we conclude 
that the differences in the number counts found between 
our field and the SDF are consistent with being due to 
cosmic variance. Finally, we note that at magnitudes 
fainter than Zr = 25 the difference in number counts be- 
comes larger again, but this is simply because our data 
are slightly shallower than the SDF. 

4.2. Comparing the LBG number density in the QSO 
field with that in the SDF 

We detect 7 objects as z ~ 6.4 galaxy candidates using 
the color criteria in § 12.71 while 9 objects are rejected 
as lower redshift objects at z < 5.5 based on a z' — Zr 
color cut (it would have been difficult to eliminate these 
objects using only i' — z'). For comparison, we apply the 
same color criteria to the SDF catalog after applying our 
magnitude limits and find only one candidate. Thus the 
number of z ~ 6.4 galaxies in the QSO field is 7 times 
larger than that the SDF field. We note that if we apply 
a more relaxed cut based on a single color of i — z > 1.9 
(the standard i-dropout technique) to the reference field 
(SDF), we still find 4 objects although 16 objects found 




Fig. 4. — Thumbnail images of the LBG candidates in the QSO 
field. The properties of the individual sources are listed in Table 

m 

in the QSO field, implying that even with the standard 
i dropout technique we find an overdensity in the QSO 
field. 

The properties of the objects identified are summarized 
in Table [TJ and we show the filter thumbnail images in 
Fig. 2J Although our data are shallower than the SDF, 
we have detected more LBGs than in the SDF. The re- 
sult suggests that the number density of z ~ 6.4 galaxies 
in this QSO field is larger than that in the SDF field. 
Such an overdensity would be consistent with the predic- 
tions from cosmological simulations (e.g., Springcl ct al. 
(2005)) that suggest that the first QSOs are situated 
in the most prominent dark matter haloes and are sur- 
rounded by a large number of fainter galaxies. Although 
the redshift discrimination offered by our photometric 
selection is not very accurate, the chance of finding a 
random projection of seven galaxies at z ~ 6 — 7 is quite 
unlikely given their overall low number density. We spec- 
ulate that we are seeing the most distant proto-cluster 
centered on a massive QSO at z = 6.43. 

4.3. Assessing the significance of the overdensity 
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Fig. 5. — The probability distribution of how many galaxies 
would be found within the field of view of Suprime-Cam using our 
selection criteria. The solid histogram is based on the simulation 
of galaxy formation at z = 7 while the dashed histogram is the 
poisson distribution assuming the average value is (7+ l)/2 = 4. 
The solid and the dashed arrows indicate the number of galaxies 
in the QSO field and the SDF comparison field, respectively. 

In this subsection, we assess how significant the over- 
density is. If we assume a Poissonian distribution, the 
probability of finding seven objects when one expects one 
(based on the SDF counts) is less than 0.01%. However 
this is not a fair statistic since we have only one com- 
parison field (the SDF) which has a similar size as the 
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TABLE 1 

The properties of the z ~ 6.4 candidate galaxies detected 
in the qso field. 



Spatial distribution of LBGs around z = 6.4 QSO 
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"2000 


<52000 
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FWHM 
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23:29:20.56 


-02:54:29.5 


24.62 


2.03 


0.44 


0.63 
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23:29:24.10 


-03:07:34.3 


24.64 


> 1.93 


0.54 


1.15 
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23:29:28.56 


-02:57:38.4 


24.69 


> 1.98 


0.52 


0.55 
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23:28:48.69 


-03:06:23.3 


24.88 


> 2.08 


0.48 


0.58 
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23:28:58.68 


-03:12:11.6 


24.98 


> 2.02 


0.41 


0.64 
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23:30:04.24 


-02:56:54.6 


25.00 


> 1.97 


0.38 


0.69 


7 


23:29:31.52 


-03:11:05.0 


25.11 


> 1.93 


0.32 


1.02 



QSO field. If we naively calculate the average number of 
counts expected from the two fields combined, (7+l)/2, 
we expect 4 counts on average. The Poissonian probabil- 
ity to find 7 objects is thus ~11%. However, the chance 
of finding one object in one field and seven in the other is 
much lower. This can be quantified by calculating the ex- 
pected number of counts that maximizes the probability 
of finding seven objects in one field and one in the other. 
Out of 10,000 Monte Carlo realizations we find 0.40% 
of such cases, i.e., the significance of this overdensity is 
99.6%. 

The above estimates are based on the oversimplifica- 
tion that galaxies are not clustered. Here we evaluate our 
results based on the cosmic variance determined from a 
cosmological simulation. The expected number of galax- 
ies is predicted from the cosmological smoothed particle 
h ydrodynamics (SPH), co ld dark matter (CDM) model 
of lChoi. J. Hi (|2009Ll20"ToT ). The simulated galaxies have 
appropriate colors and magnitudes, allowing us to make 
a direct comparison to our observation. In order to sim- 
ulate galaxies at z = 7, we have applied (i) a bright 
magnitude cut at zr > 24, and (ii) a detection complete- 
ness. Then, we count how many galaxies would be found 
in the field of view of Suprime-Cam. To evaluate the 
completeness of our detection procedure as a function 
of magnitude, we add artificial objects to the Zr image 
with all detected objects masked. The artificial objects 
are modeled by a moffat function with a half light radius 
of 0.9 arcsec in the magnitude range of 24 to 25.46. We 
used the IRAF task "noao.artdata.mkobjects" for this 
procedure. Then, we repeat the object detection and 
measure the success rate. 

After applying the completeness and magnitude cuts 
to the simulated data, we have a realistic sample that 
can be compared to our observational data. By randomly 
sampling the Suprime-cam size field of view from the sim- 
ulation, we measured the galaxy frequency distribution 
(Fig. [5]). Based on this simulation, we found that the 
probability to find one pointing with seven galaxies and 
another with one is less than 0.4%, in good agreement 
with the results from the Poissonian statistics detailed 
above. Our results suggest that it is hard to explain the 
overdensity by a chance coincidence of galaxies due to 
cosmic variance. 

4.4. Spatial distribution of the LBGs 

In this subsection, we discuss the peculiar spatial dis- 
tribution of the LBGs. Fig. [5] shows the spatial distri- 
bution of the z ~ 6.4 galaxies. The cross located at the 
center indicates the position of CFHQS J2329-0301, and 
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Fig. 6.— Spatial distribution of z > 6.4 LBGs (circles), CFHQ 
J2329-0301 at 2 = 6.43 (cross) and a newly identified QSO candi- 
date (asterisk). 

circles represent the z ~ 6.4 galaxies. It appears that 
while these galaxies are clustered around the QSO, none 
are in its direct vicinity, (the closest galaxy being already 
at a (projected) distance of ~2 Mpc from the QSO). To 
quantify the significance of this effect, we show a his- 
togram of the angular distances from the QSO to each of 
the galaxies in Fig. [7] (red solid histogram). The distri- 
bution indeed suggests that there is a deficit of galaxies 
at <2 Mpc, while typical clustered distributions are ex- 
pected to peak at the center. The blue dotted line in Fig. 
is based on the positions of randomly distributed ob- 
jects (but avoiding the masked regions just as in the real 
data). Comparing the two histograms suggests that in 
the QSO field the number of galaxies at the projected dis- 
tance from 2 to 4 (physical) Mpc/h is indeed larger than 
that expected from a random distribution. We perform 
a Kolmogorov-Smirnov test by calculating the following 
value: 



D = max I Ri — Mi 



(4) 



where Ri is the number of galaxies in the i-th bin and Mi 
is the number in the same bin as given by the randomly 
distributed galaxies. We construct a histogram of the 
D statistic using 10000 times realizations, finding that 
the distribution in the QSO field differs from the random 
distribution at a 97% confidence level. If we just compare 
the peak at 3 Mpc, it has a 1.7 a excess over the random 
distribution. 

4.5. Evidence for additional QSOs? 

Finally, we search for additional QSOs in the field by 
applying a set of color-color criteria tuned to the colors 
expected for QSOs at z ~ 6.4. Since a QSO continuum 
decreases with wavelength blue-ward of Lya, a QSO will 
have i' — z' > 1.9 and z' - z R < 0.3 (Fig. [TJ). Note that 
this is a different (bluer) cut than that used for galaxies. 
We apply these criteria to all objects detected at >3er in 
all bands and having a FWHM smaller than 0.65 arcsec 
(equal to FWHM + ctfwhm)- We found two objects (see 
Table©. One of these is CFHQS J2329-0301, showing 
that our QSO selection works. Although the redshift 
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Fig. 7. — Histogram of radial distances (in projected, physical 
Mpc) from the QSO to each of the galaxies (red line). The blue 
dashed line was derived by taking the average of 100 realizations 
of a randomly distributed population having the same number of 
galaxies. Error bars in dicate the 1 sigm a poisson error using Eqn. 
(9) and Eqn (13) from Gchrcls (1986). We find a 1.7a excess in 
the number of galaxies at a distance of ~ 3Mpc //i7o • Interestingly, 
this di stance is similar to the size of the HII region around z ~ 6 
QSO HWvithe et al.ll2005h . 



TABLE 2 

A NEW QSO CANDIDATE FOUND BASED ON COLOR-COLOR 
SELECTION. 
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4.11 
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and nature of the second object must be confirmed with 
spectroscopy, it may be another QSO associated with the 
structure of galaxies surrounding CFHQS J2329-0301. If 
so, this would add supporting evidence that this field is 
relatively overdense. The position of the QSO candidate 
is indicated by the asterisk in Fig|51 

5. DISCUSSION 

We found an overdensity of LBGs around the QSO 
with a 99.6% significance. We furthermore found evi- 
dence for a ring-shaped distribution, albeit at the <2 
sigma level. Although the physical interpretation must 
await verification using deeper data or spectroscopic 
follow-up, below we will discuss possible physical in- 
terpretations under the assumption that our results are 
significant and represent real structure surrounding the 
QSO. 

5.1. Comparison with previous work 

At lower redshift, the environments of QSOs and 
fainter A GN have been exten sively studied. At 0.05 < 
z < 0.095. iMiller et al.1 (j2003f ) found that the fraction of 
galaxies with an AGN is independent of the local galaxy 
density, in a stark contrast to both star-forming and 
passi ve galaxies that sh o w an environm e ntal d ependence 
(e.g.. IGoto et all 120031 ). ILietzen et all (|2009t ) found an 
underdensity of bright galaxies at a fe w Mpc scale from 
nearby QSOs at 0.078 < z < 0.172. iKauffmann et al.l 
(2004) found that at fixed stellar mass the number of 
galaxies that host active galactic nuclei (AGNs) with 
strong [O III] emission indicating strong AGN activity 
is twice as high in low-density regions compared to high 



density regions at 0.04 < z < 0.06. 

However the situation is different at higher redshifts. 
A Keck survey of fields centered on known z > 4 QSOs 
found excesses in the number of c o mpan ion galaxies 
(|Djorgovskilll999h . iKashikawa et al.1 (|2007j ) found that 
TBGs without Lya emission form a filamentary struc- 
ture near a QSO at z ~ 5, while Lya emitters are 
distributed a round it but avoid it within a distance of 
~ 4.5 Mpc. iMilev et al.l (|200l also found that LBGs 
are concentrated around a luminous radio galaxy at 
z « 4.1 previously found to be associated with an over- 
density of ~30 spectroscopica lly c onfirmed Lya emitters 
(iVenemans et all 120071 120021 ). In lOverzier et al.l (|2006l 
2008) it was shown that the environments of some ra- 
dio galaxies at z — 4 — 5 appear richer than the average 
field at ~ 3 — 5cr significance based on a detailed com- 
parison with GOODS. However, these studies found no 
difference between the physical properties of galaxies in 
protoclusters compared to those in the field. 

The difference between the low-z and high-z environ- 
ments of AGN perhaps stems from the different halo 
masses that host them. Due to the flux limited nature 
of observational surveys, high-z QSOs are much more 
luminous than local AGNs, and thus, presumably em- 
bedded in a more massive halo. Indeed, a generic ex- 
pectation in most models of galaxy formation is that the 
most massive density peaks in the early universe (such 
as QSOs and ma ssive galaxies) are likely t o be strongly 
clustered (|Kaiserlll984tlEfstathiou fc Reeslll988l ). In the 
hierarchical for mation and evolution scenari o of galax- 
ies and QSOs (jHaehnelt fc K auffmann 2000), luminous 
QSOs are located in rare overdense regions. This is why 
high-z QSOs ha ve been used as beacons to s earch for high 
densit y regi ons dColdwell fc La mbas 20061 ). [C room et aLl 

(2005) and iShen et all (|2007| ) found an increasing clus- 
tering of QSOs with redshift. 

Surprisingly, at z > 6 the situation appears to change 
again; luminous QSOs are n o t nece ssarily found in strong 
density peaks. iKim et aLl (|2009l ) studied the number 
densities of i-dropout objects around 5 SDSS z ^6 
QSOs using the HST/ACS. They found a n overdensity 
in two field s and underdensity in two fields. iZheng et al.1 

(2006) and Stiave llTet al.l (|2005l ) found an overdensities 
around the QSOs SDSSJ0836+0054 (z = 5.82) and SDSS 
J1030+0524 (z = 6.28). However, these overdensities 
do not appear to be as magnificent structures as were 
found at lower reds hift, or even in some random regions 
at z ~ 5 — 6 (e.g.. lOuchi et al.l[20~04bl : IKashikawa et al.l 
l2007t lOtaet al.l 120081) . 

One of the major problems of previous work at z ^6 
is the relatively small FoV of the HST/ACS. The 
200" x 200" field can only probe a region of 1 Mpc x 
1 Mpc at z ~6, and thus may easily miss any larger 
structures such as found in this work. Computer sim- 
ulations also predict that the larg est structures presen t 
at z ~ 6 span several tens of Mpcs (|Overzier et "aT] [2009') . 
while LFs of LBGs show a large fi eld-to-field variation 
(|Hu et all 12001 lOuchi et alJl2009h . These results sug- 
gest that previous non-detections of overdensities around 
QSOs may need to be re-examined using a larger field of 
view or deeper observations. The difference in the FoV 
is perhaps the primary reason why no previous work has 
detected a highly significant overdensity around any of 
the z ^6 QSOs. Our positive detection may have been 



facilitated by the large field of view that allowed us to 
investigate the structure at scale of ^3 Mpc. 
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Fig. 8. — Illustration of a massive protocluster region selected 
from a large cosmological TV-body simulation. The protocluster 
region shown is the progenitor of the most massive cluster found in 
the Millennium- 1 I Simulations (corresponding to a M ~ 10 Mq 
cluster at z = 0). Symbols show the (projected) distribution of 
dark matter halos having masses as indicated in the figure legend 
(values indicate the logarithm of the halo mass in Mq ) . The scale 
bar indicates a physical radius of 1 Mpc. 

For an illustration of what such a region at z > 6 
might look like, we show an example of a large proto- 
cluster selected fro m the Millennium-II cosmol ogical N- 
body simulations (jBovlan-Kolchin et al.l I2009I ). Figure 
[5] shows the (projected) spatial distribution of strongly 
clustered dark matter halos associated with a protoclus- 
ter at z = 6.2. This protocluster is the progenitor of the 
most massive cluster found in the Millennium-II simula- 
tions, and corresponds to a M ~ 10 15 Mq cluster when 
evolved to z — 0. Figure [3] shows, at least qualitatively, 
that we may expect significant enhancements in the den- 
sity of the galaxies that are hosted by the dark matter 
halos shown, provided that QSOs indeed trace protoclus- 
ters. A more quantitative analysis of the surface density 
of z-dropout galaxies near QSOs expected in the simula- 
tions is needed for a proper comparison. 

5.2. Size of the HII region 

Next, we compare our finding of an overdensity of 
LBGs around a QSO at a characteristic scale of 3 Mpc to 
spectroscopic measurements of HII regions around QSOs. 

It has been demonstrated through spectroscopic ob- 
servations that there exist large, ionized regions around 
luminous QSOs. These are sometimes called HII regions, 
Stromgren spheres, or highly ionized near zones. In this 
work, we will refer to them as HII regions. I n order 
to define the size of the HII region, iFan et al.l (l2006bD 
proposed a definition as a point in the spectra where 
the Lya transmission first drops to T < 0.1 for spectra 
binned in 20A pixels. The CFHQS J2329-0301 trans- 
mission d rops at T < . 1 first at 3.6 Mpc, then again at 
6.3 Mpc (|Willott et alj 120071) . Interestingly, the size of 
this spectroscopically measured HII region is comparable 



to the possible ring shape distribution of LBGs around 
CFHQS2329-0301 (Fig. Note that the size of the 

spectroscopic HII region is expected to be larger (~10 
Mpc) for the SDSS QSOs due to their higher luminosi- 
ties (they are brighter by by ~ 2 mag). In addition to 
the small FoV of HST/ACS, this may be an additional 
reason, why no significant overdensity of LBGs has been 
found around SDSS QSOs. To observe an overdensity of 
galaxies at a scale of >10 Mpc (such as expected around 
the luminous SDSS QSOs at z ~ 6), one needs multiple 
FoVs even with the Suprime-Cam. We conclude that the 
size of the HII region is consistent with the apparent lack 
of LBG candidates closest to the QSO. 

5.3. Possible Physical Mechanisms 

If our detection of the lack of LBGs near the QSO is 
real, then, what created the observed paucity of galax- 
ies within 3 Mpc from the QSO? One possibility is that 
the intense emission of ionizing radiation associated with 
QSOs ionizes the surrounding IGM and may even photo- 
evaporate the gas in neighboring dark halos before it has 
the opportunity to cool and form stars. In this scenario, 
QSOs would suppress galaxy formation in their vicin- 
ity. One would then observe a paucity of galaxies near a 
QSO despite the un derlying excess of dark matter halos. 
IShapiro et al.l (2004) presented the first theoretical simu- 
lations of the gas dynamics coupled with radiative trans- 
fer, showing that an ionizing source that emits 10 56 pho- 
tons s _1 (appropriate for a QSO) can indeed photoevap- 
orize minihalos of ~10 Mq (that would have been able to 
form small galaxies) on 100-150 Myr times cales (but see 
iBarkana fc Loebl fl999t IWvithe et al.l [2005T ) . Compared 
with our observational results of finding a "ring" -shaped 
distribution of LBGs, the central QSO may have sup- 
pressed the formation of surrounding galaxies thereby 
creating a paucity of galaxies in its vicinity but still hav- 
ing an overdense region beyond the inner, ionized region. 
This may explain the observed results, at least qualita- 
tively. 

However, iKashikawa et al.l (|2007[ ) quantitatively esti- 
mated that such UV radiation from QSO can suppress 
galaxy formation or not. According to their simulation, 
QSO UV radiation can suppress star formation in a halo 
with M vir < 1O 9 M , while a halo with M mr > 10 n M Q is 
unaffected. Although mass estimate of our z ~ 6.4 LBG 
are uncertain because we do not have deep near-infrared 
data, considering bright magnitude, they are likely to be 
massive galaxies with M vir - 10 11_12 M Q . If so, QSO 
UV radiation is not strong enough to create the central 
deficit we observed. 

Can we find other theoretical predictions that 
could explain the obs ervations?. The simulation of 
IKashikawa et al.l (|2007l ) assumed star formation in spher- 
ical dark matter halos. They found that the impact of 
photoionization is greater if that star formation is tak- 
ing place in a disk or in substructures (which is perhaps 
likely to be the case). If the QSO formed at a much 
earlier time than the surrounding galaxies, this scenario 
might be able to suppress the formation of galaxy seeds 
at the time when their masses were still sufficiently small. 

Another plausible scenario is that the IGM surround- 
ing the QSO may have been ionized by galaxies long be- 
fore the QSO turned on. It has been suggested that be- 
cause luminous QSOs are expected to form in rare over- 
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dense regions, the surro unding IGM ha d already been 
pre-ionized by g alaxies (lYu fc Lul 120051) . According to 
the simulation of lBolton fc: Haehneltl (|2007T ). the neutral 
hydrogen fraction, fHi, near a QSO is estimated to be 
fm < 0.3. Therefore, QSO radiation is emitted into a 
substantially pre-ionized IGM. In this case, galaxies that 
formed before the QSO started emitting its strong UV 
radiation may still be present in the vicinity of the QSO. 
Because such galaxies are likely to have ceased their star 
formation, they would not be detected by our LBG tech- 
nique. 

In an overdense region, the QSO host galaxy is likely 
to be the mo st massive galax y that formed the earliest. 
According to lYu fc Lul (|2005| ). QSO host galaxies expe- 
riencing rapid star formation at a rate of ~3000M© yr _1 
combined with the radiation field emitted by the QSO 
itself can produce enough photons to ionize a large HII 
region. In order to check whether we see any evidence for 
the presence of such a massive host galaxy, we have per- 
formed a careful PSF subtraction on our images. The 
residuals indicate the detection of a large host galaxy 
by SED fitting (i? e =llkpc, < 10 10 M^) associa ted with 
the QSO CFHQS J2329-0301 (IGoto et al.ll2009l) with ev- 
idence of extensive star formation based on its extended 
rest-frame UV flux. The presence of a massive host 
galaxy may thus support a scenario in which the QSO 
and its host galaxy evolve together, suppressing galaxy 
formation in their vicinity. 

Alternatively , it has also been suggested by 
iDiikstra et all (|2008l ) that in order to facilitate the 
formation of a supermassive black hole by z ~ 6 in the 
first place, it may be required to have a rare pair of 
dark matter halos in which the intense UV radiation 
from one halo prevents fragmentation of the other, so 
that the gas collapses directly into a supermassive black 
hole, and explaining the lack of galaxies forming in its 
vicinity as observed in our field. Consistent with this 
scenario, it is interesting to note that we found a second 
QSO candidate located near the eastern edge of the 
field (see §4.4). We speculate that this could be a QSO 
associated with the adjacent (second) halo that ionized 
the original halo. 

Last, it has been suggested that az~6 QSO is likely 
to have experienced multiple mergers in order for its 
black hole to grow to a mass of ~ 1Q 9 M(^ . An example o f 
this scenario is given by the simulation of lLi et ail (|2007f ). 
which predicts that the host galaxy of a z ~ 6 QSO may 
have experienced seven mergers with mass ratios of 4:1 
or greater (see their Fig. 3). The QSO in our field may 
thus have merged with all the galaxies in its direct vicin- 
ity, explaining the peculiar spatial distribution that we 
find. This is also consistent with the disco very of the 
large host galaxy associated with the QSO ()Goto et al.l 
[20091) . 

In summary, although our detection of substructure 
in the LBG distribution near the QSO is rather weak 
(< 2a), studies of the spatial distribution of LBGs 
around QSOs at z ~ 6 are very important for testing 
numerous of the interesting physical scenarios related to 
the co-evolution of QSOs and galaxies as discussed above. 
Therefore, it is important that our conclusions are veri- 
fied using deeper, multi- wavelength imaging and spectro- 
scopic observations that may detect additional galaxies 



missed by our current selection, such as galaxies that are 
below our (UV) detection limits, star-forming galaxies 
that are heavily obscured, or galaxies in which star for- 
mation has ceased. Also, it will be important to extend 
the analysis performed here to other fields to obtain good 
statistics on the typical structures associated with QSOs 
at this extreme redshift. 

6. CONCLUSIONS 

Taking advantage of the large field of view (34'x27') 
and the new red-sensitive fully depleted CCDs recently 
installed on the Subaru/Suprime-Cam, we have investi- 
gated the large scale environment around the most dis- 
tant QSO studied to date. 

Our findings are as follows. The number of candidate 
LBGs at z ~ 6.4 is 7 times larger than that in a com- 
parison field (the Subaru Deep Field), suggesting that 
the QSO field hosts an overdense region. We estimate 
that the probability that this overdensity is a chance co- 
incidence is less than 0.4% based on either using simple 
Poissonian statistics, or on a comparison with a cosmo- 
logical SPH simulation of z ~ 7 galaxies. 

We find evidence for a non-uniform distribution of the 
LBGs in a "ring-like" shape surrounding the QSO at a 
radius of ^3 physical Mpc, i.e., galaxies are overdense 
around the QSO, but at the same time they avoid the 
very center near the QSO (see Figs. [6] and [7]). A KS- 
test shows that the radial distribution of LBGs in Fig. 
[7] is different from random at a 98% significance level. 
The distance of 3 Mpc is compara ble to the size of th e 
HII region around QSOs at z ~ 6 (jWvithe et al.l 120051) . 
Possible physical explanations of such a central deficit of 
galaxies include the suppression of galaxy formation due 
to the strong UV radiation field of the QSO. However, be- 
cause the significance of our detection of a non-uniform 
distribution of LBGs around the QSO is low, it is im- 
portant to verify these results with deeper imaging data 
and/or spectroscopy. 

Our findings show that QSOs at z ~ 6 may indeed be 
embedded in the densest regions of the early universe, 
provided that they are observed on significantly larger 
scales compared to previous studies that used the rela- 
tively small field of view provided by HST. 
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